function [freqs, resp] = etffft(u, y, T0)

%untested. supposed to work well with a chirp signal without stochastic
%components. returns G(exp(i*freq*T0) = G(exp(i*2*pi*(1:(N/2))/N), where G
%is the z-space transfer function.

N = length(u);
if(length(y) ~= N) 
    error('vector lengths do not match');
end

freqs = 2*pi*(1:(N/2))/(N*T0);
resp = fft(y)./fft(u);

end